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Shower spectrum and detection of tau from Earth-skimming neutrinos 
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A realistic Monte-Carlo simulation code is developed to study the shower spectrum initiated by decay of tau 
leptons, which come from charged current interaction inside the Earth. AGN and GZK neutrino fluxes are used 
to simulate the shower spectrum. 
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At neutrino energy above lO^^eV, detecting 
neutrinos skimming through the Earth provides 
an unique chance to detect those high energy neu- 
trinos and confirmation of ~^ Vt oscillation. To 
evaluate the detectable event rate and optimize 
the detection efficiency, a realistic Monte-Carlo 
simulation code is developed. 

The Earth is modeled as a spherical sphere of 
several layers of materials and density. The v — N 
cross-section is modeled by CTEQ-6 model J]. 
For this study, the energy loss of r is modeled by 
continous energy loss, i.e., 

-dE/{pdx) = a + /3E = f3E, 

where p is density of material, dx is distance trav- 
eled by T, and a and /3 is the coefficients of energy 
loss from all processes PP|. The r survival proba- 
bility P is modeled by 

dP/P = ~dx/{KE) = dE/{K(3pE'^), 

where kE is decay length of r at energy E in PeV 
and K = 49.9 m/PeV. 

T decay is simulated with TAUOLA package 
with fully polarized ts. When ts decay inside 
the Earth, the i^tS are re-propagated inside the 
Earth and could generate new ts again. When a 
r decays outside the Earth, it initiate air shower, 
except when r decay to muon. The shower energy 
is sum of the energy from electrons and hadrons. 
The shower simulation is done by a modified ver- 
sion of CORSIKA 12]. 
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To compare with previous results from analyt- 
ical method UrS are injected to a fix width 
of 100km of rock. Shown in the Fig. 1, the r 
spectrum from both methods agrees quite well. 
This test conclude that the algorithm used in this 
Monte-Carlo simulation is correct. We are updat- 
ing energy loss with Stochastic processes and the 
complete study will be published soon. 
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Figure 1. Spectrum of r from i>t- pass through 100 
km of rock from AGN (left figure) O and GZK 
(right figure) 0. The blue markers come from 
analytical calculation PP, while the red markers 
come from Monte-Carlo simulation in this study. 
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